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We investigate both analytically and numerically the ground-state and thermodynamic 



properties of the quantum mixed spin-1/2-1/2-1-1 chain described by the Hamiltonian 
H = J2t=i (Ji $a-3 ■ SU-2 + J2 su-2 ■ o4t-i + J3 Su-1 ■ <S 4 £ + J 2 Su ■ su+i) , where two £=1/2 
*vj ' spins (s4£_3 and 84^—2) and two S—l spins (Sa-i and So) are arranged alternatively. In 

several limiting cases of J%, J2, and J3 we apply the Wigner-Eckart theorem and carry out a 
perturbation calculation to examine the behavior of the massless lines where the energy gap 



vanishes. Performing a quantum Monte Carlo calculation without global flips at a sufficiently 
low temperature for the case where J\ = J3 = 1.0 and J2 > 0, we find that the ground state of 
the present system in this case undergoes a second-order phase transition accompanying the 
vanishing of the energy gap at J 2 = 'he with J 2c = 0.77 ± 0.01. We also find that the ground 
states for both J 2 < J 2c and J 2 > J 2c can be understood by means of the valence-bond-solid 
picture. A quantum Monte Carlo calculation which takes the global flips along the Trotter 
*^J ■ direction into account is carried out to elucidate the temperature dependences of the specific 

heat and the magnetic susceptibility. In particular, it is found that the susceptibility per unit 
cell for Ji = 0.77 with Ji = J3 — 1.0 takes a finite value at absolute zero temperature and that 
the specific heat per unit cell versus temperature curve for J 2 = 5.0 with Ji = J3 = 1.0 has a 
double peak. 
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§1. Introduction 

A quantum antiferromagnetic chain has been the subject of numerous experimental as well as 
theoretical studies over a long period, in particular since Haldane's prediction!!!' implying that 
the chain with an integer spin is quite different from the chain with a half-integer spin in the 
ground-state and low-lying-excited-state properties. Almost all results of these studies support the 
prediction, and now it is widely agreed that for the case of isotropic nearest-neighbor exchange 
interactions, the former chain has massive excitations with a exponential decay of the two-spin 
correlation functions, while the latter chain has massless excitations with a power-law decay of the 
two-spin correlation functions Q> 

Under these circumstances, a mixed spin chain consisting of two kinds of quantum spins which 
are arranged periodically has been of considerable interest in recent years. On the experimental 
side, several real materials of this type have been synthesized, and their magnetic properties have 
been observed!!) A typical example of these materials is CuNi(pba)(H20)3-2H20 with pba=l,3- 
propylenebis(oxamato),o) where an S = l/2 spin of the Cu 2+ ion and an5 = l spin of the Ni 2+ ion 
make an alternating array and each nearest-neighbor pair of the spins couples by the isotropic, anti- 
ferromagnetic exchange interaction. On the theoretical side, on the other hand, a variety of models 
describing the quantum mixed spin chain have been investigatedBB0&B , ll3'llJ'll3t3EJ'Ej'llB' Inte- 
grable models have been constructed and discussed by several authors .ErBLrcP Although these inte- 
grable models are specific ones in which rather complicated interactions are assumed, exact solutions 
for the models are helpful for us to understand the essential consequences of the quantum mixed 
spin chain. Motivated by the above experimental observations, several authorsoE^EJ'E3ESE§E!P 
have studied both analytically and numerically a simpler case of the 5 = 1/2 and 5 = 1 alternating 
spin chain with antiferromagnetic nearest-neighbor exchange interactions. A characteristic feature 
of this chain in the case of isotropic or Ising-type interactions is that it is a quantum ferrimagnet 
and its ground state is magnetic. 

Generally speaking, quantum fluctuations play a more crucial role in the nonmagnetic ground 
state than in the magnetic one. As an example of the quantum mixed spin chains which have 
isotropic nearest-neighbor exchange interactions only and also have a nonmagnetic ground state, 
we consider in this paper a one-dimensional antiferromagnet where two 5= 1/2 and two 5 = 1 spins 
are arranged alternatively. For the sake of simplicity we call this system the spin-1/2-1/2-1-1 
chain. The Hamiltonian is given by 

N/A ( \ 

H = 2_^ ( Jl S*4f-3 - •?4£-2 + Jl Sii-2 ' BU-\ + ^3 SW-1 ' &U + J2 &U ' •?4£+l J , (1-1) 

where S£ with £ = 1 and 2 (mod 4) and Si with £ = 3 and 4 (mod 4) are, respectively, the 5 = 1/2 
and 5 = 1 spin operators at the £th site; J\, J2, and J3 are, respectively, the interaction constants 



between the nearest-neighbor pair of 5 = 1/2 spins, that of 5 = 1/2 and 5 = 1 spins, and that of 
5 = 1 spins; N, being assumed to be a multiple of four, is the total number of spins. We impose 
periodic boundary conditions (sjv+i = si). 

It is naturally anticipated that the present mixed spin system has different ground states depend- 
ing upon the signs of J\, J2, and J3. As we will discuss later in more detail (see Table I below), 
the value <S t ot of the total spin 

N/4 
5tot = X, 1^41-3 + S4£_2 + Su-l + Suj (1-2) 

in the ground state can be determined rigorously by applying the Lieb-Mattis theorem.tiP Except in 
the region where J\ < and J3 < 0, the ground state belongs to the 5tot = subspace, and therefore it 
is nonmagnetic. Furthermore, in the region where J\ < 0, J2 > 0, and J3 < and in the region where 
J\ <0, J2 <0, and J3 <0, the ground states belong, respectively, to the <S to t =^/4 and 5 to t =3iV/4 
subspaces. The former is ferrimagnetic, while the latter is ferromagnetic. For comparison we depict 
in Fig. 1 the spin array with the lowest energy of the corresponding classical system in which si and 
Si are replaced, respectively, by the classical spin vectors whose magnitudes are 1/2 and 1. From 
this figure and Table I we find that there is a close analogy between the classical lowest-energy state 
and the quantum ground state. It is interesting to see that the periods concerning the translational 
symmetry of the classical lowest-energy states are, respectively, four and eight in units of the lattice 
spacing in the region where JiJ^>0 and in the region where JiJs<0. 

We explore here both analytically and numerically the ground-state and thermodynamic prop- 
erties of the present system]!^ In several limiting cases of the interaction constants we apply the 
Wigner-Eckart theoremlij' and carry out a perturbation calculation to examine the behavior of the 
massless lines in the J\ versus J% plane with Ji>0 and J3>0. Confining ourselves to the case 
where J\ = J3 = 1.0 and J2>0, we perform numerical calculations. In these calculations we mainly 
employ a quantum Monte Carlo (QMC) method, and do an exact-diagonalization calculation for 
iV = 8 and 16 to check the QMC results. To determine as accurately as possible the massless point 
J2 C of J2 in the above case is one of our main purposes. We discuss how the ground states for both 
J2 < <^2c and J2 > J2C are represented in the frame of the valence-bond-solid (VBS) picture proposed 
by Affleck, Kennedy, Lieb, and TasakiHD We also aim at clarifying the temperature dependences 
of the specific heat and the magnetic susceptibility for a few values of J2 in this case. Very recently, 
Fukui and KawakamiliiP have also discussed this system. They have mapped this system to the 
non-linear sigma model and have shown that in its ground state a second-order phase transition 
which accompanies the vanishing of the energy gap may occur. 

The remainder of the paper is organized as follows. The next section (§2) and the following 
one (§3) are devoted to analytical treatments and numerical calculations, respectively. Finally, the 
main results are summarized and further discussed in §4. 



§2. Analytical Treatments 

We start with applying the Lieb-Mattis theoremliiP to determine the value <Stot of the total spin 
<Stot i n the ground state for arbitrary values of J\, J2, and J3. In the present case the theorem 
implies that, if we can divide the lattice into the A and B sublattices in such a way that all the 
intrasublattice interaction constants are not positive and all the intersublattice interaction constants 
are not negative, then «S tot is given by 



-Hot 



'max, A <->max,i} 



(2.1) 



where 5 max ,A and S ma x,B are, respectively, the maximum values of the total spins in the A and 
B sublattices. (Note that the constant g 2 in ref. 17 is equal to zero because the Hamiltonian Ji 
consists of the nearest-neighbor interactions only.) Let us consider two examples. For the region 
where J\ >0, J2 >0, and J3 >0, we choose the sites £a belonging to the A sublattice to be I a = 1 
(mod 2) and the sites £3 belonging to the B sublattice to be £3 = % (mod 2). Then, all of Ji, J2, 
and J3 are the intersublattice interaction constants, which satisfy the above requirement, and thus 
we obtain 5 tot = since <S max ,A = 5 max ,,B = 3./V/8. For the region where J\ > 0, J2 < 0, and J3 > 0, on 
the other hand, we choose I a to be £a = 1 and 4 (mod 4) and £3 to be £3 = 2 and 3 (mod 4). Then, 
J2 is the intrasublattice interaction constant and the remaining J\ and J3 are the intersublattice 
interaction constants, which again satisfy the above requirement. Thus, we obtain again 5t t = 
since 5 maXj ^ = 5 maXj B = 3iV/8. In a similar way, we can determine <St t f° r the remaining regions of 
Ji, J2, and J3. The results are listed in Table I together with £a, £b, Sm&x,A, and S maXj B for each 
region. It should be noted that the A and B sublattices correspond, respectively, to the 'up' and 
'down' sublattices in the lowest-energy classical-spin arrays shown in Fig. 1. 

We now turn to an application of the Wigner-Eckart theoremllSP to the following three limiting 
cases: (a) the case where | J2 1 S> | J\ \ , | J3 1 with J2 > 0, (b) the case where | J2 1 S> | J\ \ , | J3 1 with J2 < 0, 
and (c) the case where |Ji|3>|«/2|) |«^3| with J\ <0. In the limiting case (a), a pair of the spins s^_2 
and S4£_i can be replaced by an effective S' = l/2 spin •s' 4 < | f l 24 ^_ 1 . This replacement means that we 
take into account only the doublet state of the two-spin system consisting of su-2 and Su-i and 
neglect the quartet state of the system. Similarly, a pair of the spins Su and su+i can be replaced 
by an effective 5 = 1/2 spin s£f Ae+1 Then, the Wigner-Eckart theorem implies that both s*4^_2 and 
S41-1 are proportional to s^_ 24 £_i and both S41 and s*4£ + i are proportional to s^ 4£+1 : 

su-2 = ol s 4 e £ _2,4f-i > (2.2a) 

S41-1 = P s£i-2,4:e-i i (2.2b) 

Su = (3' 541,41+1 ) (2.2c) 

sa+i = a' slf M+1 , (2. 2d) 



Table I. Values of £a, Cb, S max ,A, <S max ,s, and <S to t for each region of Ji, J 2 , and J 3 . 



Region of Ji , J 2 , and J 3 


£4 


£b 


'-'max, A 


«->max,B 


Stot 


Ji> 0, J 2 >0, J 3 >0 


1 (mod 2) 


2 (mod 2) 


W/8 


3iV/8 





Ji> 0, J 2 <0, J 3 >0 


1, 4 (mod 4) 


2, 4 (mod 4) 


3iV/8 


3iV/8 





Ji< 0, J 2 <0, J 3 >0 


1, 2, 3, 8 (mod 8) 


4, 5, 6, 7 (mod 8) 


3AT/8 


3iV/8 





Ji< 0, J 2 >0, J 3 >0 


1, 2, 4, 7 (mod 8) 


3, 5, 6, 8 (mod 8) 


37V/8 


3iV/8 





Ji> 0, J 2 >0, J 3 <0 


1, 3, 4, 6 (mod 8) 


2, 5, 7, 8 (mod 8) 


37V/8 


3AT/8 





Ji> 0, J 2 <0, J 3 <0 


1, 6, 7, 8 (mod 8) 


2, 3, 4, 5 (mod 8) 


3AT/8 


3N/8 





Ji< 0, J 2 <0, J 3 <0 


all 


none 


3AT/4 





37V/4 


Ji< 0, J 2 >0, J 3 <0 


1, 2 (mod 4) 


3, 4 (mod 4) 


N/4 


N/2 


N/4 



a, (3, (3' , and a' being constants. The values of a, (3, (3' , and a' can be obtained, respectively, by 
calculating all matrix elements of both sides of eq. (2.2a), eq. (2.2b), eq. (2.2c), and eq. (2. 2d) in 
the doublet subspace. The results are a = a' = —1/3 and (3 = (3' = 4/3. Substituting eqs. (2.2a)- 
(2. 2d) with these constants into eq. (1.1), we obtain, apart from the constant term coming from the 
J2-terms, 

T-l — ST I ! T 5* cff ? cff _L 16 T ? cff ? cff I <0 f\ 

rl ( a ) ~~ Z^ I q l s 4£-4,M-3 ' S U-2,U-1 "+" "77" J 3 s 4^-2,4^-l ' s 4£,4£+l I K*-'*) 

with So^^^JVAf+i — ^jVI- Thus, we can conclude that in the limiting case (a), our system described 
by the Hamiltonian TL of eq. (1.1) is equivalent to the 5 = 1/2 bond-alternating chain described by 
the Hamiltonian Hfa) of eq. (2.3), as far as the ground state and the sufficiently low-energy excited 
states are concerned. 

In the limiting case (b), we replace a pair of the spins su-2 and <S4£_i {Su and s*4^+i) by an 
effective S = 3/2 spin S 4 ( |l 24 ^_ 1 (S 4 C ^ +1 ), neglecting the doublet state of the two-spin system 
consisting of s*4^_2 and S^-i (S41 and s*4^ + i). Applying the Wigner-Eckart theorem, we obtain 

1 ""eft 
He-2 = ^A£-2,U-1 1 (2.4a) 

-> 2 -> 

SaI-1 = o S 4^-2,4£-l > ( 2 - 4b ) 

-> 2 -> 

SaI = o S 4^,4£+l ) ( 2 - 4c ) 

su+i = -^u,u+i ■ (2.4d) 




Thus, we can show that in the limiting case (b), our system is equivalent to the 5 = 3/2 bond- 
alternating chain described by the following Hamiltonian: 

"^( b ) = Z^ ( q^ 1 ^4f-4,4l-3 ' ^4£-2,4£-l + g^3 ^4^-2,4^-1 " ^U,U+1 J (2-5) 

with Sq 5 ^ = SjY ff Ar+1 = SjY ff 1 , as far as the ground state and the sufficiently low-energy excited states 
are concerned. Furthermore, in the limiting case (c), we replace a pair of the spins s*4^_3 and s*4£_2 
by an effective 5=1 spin Sf^_ 3Ae _ 2 , neglecting the singlet state of the two-spin system consisting 
of S4^_3 and su-2- The Wigner-Eckart theorem gives 

1 



feflf 
1 



2^U-3,U-2 > (2.6a) 



S41-2 — 2^41-3,4^-2 • (2.6b) 

Thus, our system in the limiting case (c) is equivalent to the 5 = 1 chain described by the following 
Hamiltonian: 

N 'U l -off - 1 - -eff 

^(c) = Z^ I o^ 2 ^4£-3,4£-2 ' Su-1 + ^3 ^-l • Su + - J2 fi« • S^ +1>u+2 ) (2.7) 

with <Sjy^. 1 N+2 = 5 1 c f , as far as the ground state and the sufficiently low-energy excited states are 
concerned. 

From the above arguments for the first two limiting cases (a) and (b), we can discuss when the 
system has massless excitations in these cases. In the case (a), it is massless only when 

Ji = 16 J 3 > , (2.8) 

since the antiferromagnetic 5 = 1/2 bond-alternating chain becomes massless only in the case of no 
bond-alternation.ElP In the case (b), on the other hand, the system is massless whenEBEj^ 

Jj = 4J 3 > and Ji = 4J 3 ^^ > . (2.9) 

1 =F(5 



The value of 5 has been deter mined" 123 recently to be 5 = .43 ± 0.01, which is the value given in 
ref. 25. 

We also investigate the possibility of massless lines in the limiting cases: (d) the case where 
J3 ^>\Ji I, [J2 1 an d ( e ) the case where J\^>\J2\ and Ji3> J3>0. For this purpose we perform 
perturbation calculations. For the case (d), we first observe that the system is massless at the point 
Ji = J2 = 0.0 when J3 > 0. This is because at this point the system consists of iV/4 independent 
pairs of two 5 = 1 spins, each of which has a singlet state as the lowest-energy state, plus N/2 free 
5 = 1/2 spins, and therefore the ground state of the system is 2 N ' 2 -fold degenerate. In order to 



explore the neighborhood of the above point for finite Ja(> 0), we consider the four-spin system 
described by the Hamiltonian, 

h e = J2 S*4£-2 • S4J-I + <^3 Su-1 ' Su + Jl Su ■ «Wl • (2-10) 



Diagonalizing analytically this Hamiltonian, we find that the ground and first-excited states of the 
four-spin system are singlet and triplet states, respectively. The energy difference J between them 
is given by 



j=m*\-t\+*\-t\ -^m +o m , (2.1D 




while the energies of the second- and higher-excited states measured from the ground-state energy 
are finite at J2 = 0.0. Thus, taking only the ground and first-excited states into consideration, we 
regard the four-spin system as the system described by JBi ■ fg, where both &i and fg are S = 1/2 
operators. Then, as far as the ground state and the sufficiently low-energy excited states are 
concerned, we map our original system described by the Hamiltonian TC of eq. (1.1) to the following 
effective Hamiltonian: 

N/i ( . _ . \ 

HeS = 53 ( J°t ' ^ + *^ 1 ^ ' ° e + 1 + ^ £ ' ° l + l + ^ 2 ^ ' ^+1 J (2-12) 

with Bn , 1 = B\ and tn ,-> = T\. The values of J\ and J2 are determined in such a way that the 

4 * „ 

matrix elements of the two operators J\Tg ■ ae+i+J2&e • (Ti+x+J-zTg ■ f^+i and J\Su + i ■ 841+2 coincide. 

The mapping is correct to the order of {J2/J3) 3 , and J\ and J2 are given by 

^i=^|l-^(|) 2 -^(^) +0|m ]}, (2.13a) 

j2=ji {Kt) 3+o "^ 1 ]y (2 - i3b) 

Thus, to the order of (^/^s) 3 ; the effective Hamiltonian 7i e s describes the S = 1/2 chain which 
has the bond-alternating nearest-neighbor interactions with the interaction constants J and J\ 
and also the uniform next-nearest-neighbor interaction with the interaction constant J2- Thus, 
remembering the fact that the antiferromagnetic 5 = 1/2 chain with the uniform nearest-neighbor 
and next-nearest-neighbor interactions where the ratio of the latter interaction to the former one 
is less than 0.24 is massless,cj) we may conclude, from eqs. (2.11), (2.13a), and (2.13b), that in 




the limiting case (d), the present quantum mixed spin chain described by the Hamiltonian fi of 
eq. (1.1) is massless when 

-^(!)+<f)^(tf+#f IJ - <2 ' 14) 

For the case (e) we perform a similar perturbation calculation in the following way. We first 
diagonalize analytically the Hamiltonian which describes another four-spin system, 

h £ = J2 Su-i ■ S41-3 + J\ S*4£_3 • S*4^_2 + J2 3u-2 • SW-1 , (2-15) 



to obtain the energy difference J between the singlet ground and triplet first-excited states and 
that J' between the singlet ground and quintet second-excited state to be 

■^•M^l +OI \~~r\ \), ( 2 - 16a ) 




J, = J iAi) +0 \\i))y (2 - 16b) 

Adding the fact that the energy difference between these lowest states and higher-excited states 
are finite at | J2I/J1 =0.0, we can regard, to the order of (J2/J1) 2 , the four-spin system as a system 
described by Jag ■ ft, where Bg and T£ are S = 1 operators. Then, as far as the ground state and 
the sufficiently low-energy excited states are concerned, we map the present system described by 
the Hamiltonian TL of eq. (1.1) to the effective Hamiltonian of eq. (2.12) with the 5 = 1 operators 5g 
and fg. The values of J\ and J2 can be determined similarly to the above perturbation calculation. 
The results are given by 

j -4"Kt) 2+o ((tf)}' (2 - l7a) 

J2=J3 °((f) )' (2 - 17b) 

which show that, to the order of (J2/J1) 2 , the effective Hamiltonian 7i c s describes the S = l chain 
which has the bond-alternating nearest-neighbor interactions with the interaction constants J and 
J\. Thus, it may be concluded, from eqs. (2.16a) and (2.17a), that in the limiting case (e), the 
present system is massless when (Ji/2)(J 2 /Ji) 2 /j 3 jl - (1/2)(J 2 /Ji) 2 | = (l ± 5) / (l =f 6) with 
5~l/4,E3'E3l 2 3Ej'Ey) or equivalently, when 



The above results given by eqs. (2.8), (2.9), (2.14), and (2.18) with J3 fixed at a positive and 
finite value describe, respectively, the asymptotic behavior in the limiting cases (a), (b), (d), and 
(e) of the massless lines on the J\ versus J2 plane. These are summarized in Fig. 2, where J3 is 
chosen to be J3 = 1.0. Our preliminary results indicate that the asymptotic massless lines in Fig. 2 
connect to form four massless lines. Leaving the detailed presentation of these massless lines for 
further publication,ElP we determine in the next section the massless point of J2 for the case where 
Ji = J-^ = 1.0 and J2>0, as mentioned in §1. 

§3. Numerical calculations 

3. 1 QMC method 

As has been mentioned in §1, we mainly employ a QMC method in the numerical calculation. 
Before discussing the numerical results, we explain the QMC procedure for the present spin-1/2- 
1/2-1-1 chain. By the use of the Suzuki- Trotter decomposition^ 2 ]) of checker-board type,E2P the 
partition function Z = Tv exp(— (3TC) , where (5= (k^T)^ 1 with the Boltzmann constant k-Q and the 
absolute temperature T, is approximated as 



Z~Tr 



N/4 ( \ 1 

II \ ex p(-/ ?/l 4£-3,4£-2/nJ exp(-/3/i 4 ^_ M ^/nJ \ 



N/4 ( 



Y[ I exp(-/?/i 4<? _ 2 ,4<?-i/™) exp(-/?/i 4<?)4m /n 



e=i (. 



(3.1) 



Here, 



hu-3,4£-2 = J\ 84£-3 • S*4^-2 , 

hu-i,u = J3 Su-i ■ Su , 

h4e-2,u-i = J2 S41-2 • Su-x , 

h4i,4i+\ = J2 Su ■ Su+1 



(3.2a) 
(3.2b) 
(3.2c) 
(3.2d) 



are the local Hamiltonians, and n is the Trotter number. Choosing the z-axis as the quantization 
axis, we introduce the local Boltzmann factors defined by 



(2r-l,2r) 
PU-3,4£-2 — 


/ (2r-l) (2r 
\ 6 4£-3 b 4i- 


-1) 
-2 


exp(-/3/i 4 £_ 3i4 £_2/ 


'n\ 


,(*) J 2r ) \ 
b 4l-Z b 4l-2 J > 


(3.3a) 


(2r-l,2r) 


_ 1 «(2r-l) <y(2r-l) 
" V °U-l °4£ 


exp(-(3hu~i,u/n) 


«(2r) o(2r)\ 


(3.3b) 


(2r,2r+l) __ - 


1 (2r) q (2r) 

b 4l-2°4l-\ 


exp 


■(■ 


-/3hu-2,4£-i/ n ) 


b 4t 


■+l)„(2r+l)\ 
-2 °U-1 J ' 


(3.3c) 



(2r,2r+l) _ / q (2r) (2r) 
P4£,4£+l ~ \ ^ S 4£+l 



exp 



(-Ph. 



UM+1 



jn 



<j(2r+l) (2r+l) 
°U b U+l 



(3.3d) 



with 



>') >') 



s 4£-3 S 4£-2/' 



o(r') o(r') 

J u-i D 4e 



>') q (r>) 



»4£_2"/«-l/' an °- 

i-2, of S4|_i an< 
Here, r' is a label along the Trotter direction, taking the values, 1, 2 
take the values, ±|; <% 
eq. (3.1) is rewritten as 



(r') tr') \ 

^le s 4£+i); which denote, respectively, the spin 



states consisting of s^_3 and s^_2, of S^-i and S^, of s*4£_2 and S^-i, and of fi^ and su+i- 

o . M j (r') 

*, Z7Z, ^4^_3 anQ ^4^_2 

take the values, ±^; 5^_ 1 and o^ take the values, 0, ±1. Using these local Boltzmann factors, 



Z 



E 



E 



E 



__ n (N/4/ v 

Z^ 11)111 #41-3,4^-2 A«-l,4£ 



{ S «-3H S 4(-2}{*4(-lH S 4( } I 



JV/4 



n| (2r,2r+l) (2r,2r+l) 
I ^4£-2,4^-l A«,4£+l 



(3.4) 



where /%-2,4<Li — /°L-2,-«-i anci Pu,u+i = pJ'u+v Tlie right-hand side of eq. (3.4) can be 
regarded as the partition function of a two-dimensional Ising system with the Ising variables, 
{ s 4£-3}' { s u-2^ {.^u-iii an d {$ai }' wn i cn has four-body interactions corresponding to the local 
Boltzmann factors. A graphical representation of this two-dimensional Ising system is presented in 
Fig. 3. It is noted that in this figure we have four kinds of plaquettes, corresponding to four kinds 
of the local Boltzmann factors. 

We perform a QMC calculation on the basis of the two-dimensional Ising system discussed above. 
To update the spin configuration, we carry out the following local flips, keeping the total magneti- 
zation of the system constant :tS 



,(2r) 



,(2r) 



S 4f-3 ~~ 2 ' S 4£-2 



1 (2r+l) 
2' S 4£-3 



£ (2r+l) 

2' S 4^-2 



,(2r) 

'4^-3 



X, S 



(2r) 



' °4£-2 



-, S 



(2r+l) 



' J 4^-3 



I ( 2r +l) _ 1 1 

2 ' 



2' S 4^-2 



(3.5a) 



,(2r-l) 



5 



(2r-l) (2r) 



4£-l ' S 4^-2 ~~ 2' 4^-1 



5 



(2r) 



,(2r-l) 



1 c(2r-l) + i g 
; "4£_i i^ r, 6 



(2r) 

4^-2 



-- 5 (2r) +ll 
2 > D 4£-l + f ' 



(3.5b) 



q (2r) q (2r) q (2r+l) <j(2r+l) 



?(2r) _l -, c (2r) _ ! Q (2r+1) _,_ 1 Q (2r+1) _ j_ 
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(3.5c) 



Q (2r) _ 1 Q (2r) 
°4£-l — x ' °4^ 



1 Q (2r+1) _ 1 c(2r+l) _ 1 

_i ) d ai-i — - 1 ' °4^ — L 



10 



o(2r) _ , o(2r) _ 1 c(2r+l) _ - „(2r+l) _ , I /o y\ 

°4£— 1 — x > *- , 4£ — -*-> °4£— 1 — x ' °4£ — I ' I^O.-JUJ 



q (2r-l) (2r-l) _ 1 Q (2r) (2r) 

<->4£ ' s 4£+l ~~ 2 ' 4 ^ ' S 4£+l 



? (2r-l) r „(2r-l) _ , 1 c (2r) _ -, „(2r) 

2' 



9 (2r-lj _ 1 (2r-l) _ , i- ^rj _ (2r; _ - I ,, _ , 



*^4£-i (^' = 2r — 1, 2r, 2r+l) in eqs. (3.5b) and (3.5c) being equal to — 1 or 0, and S^ e (r' = 2r— 1, 
2r, 2r+l) in eqs. (3.5c) and (3.5e) is equal to or 1. 
We also carry out the global flips along the Trotter direction, which are given byt^ 



_(1) _ I ,(2) _ I ,(2n) _ I I 

»// ~ 2 ,S t , - 2 , ,V - 2 j 

= -~, sf = -i, • • • , 4? n) = -~ J (£' = 4£ - 3, 4£ - 2) , (3.6a) 

J c(!) c( 2 ) c(2n)\ 

a {' ' D fJ ' ' ' ' ' °£' f 

*± 1 4 1} - 1, 4 2) - !> • • • > 4 2n) - 1 } (£' = U-1, U), (3.6b) 

4' ( r ' = 1) 2, •••, 2n) in eq. (3.6b) being equal to or 1. These global flips do change the 
magnetization in contrast to the above local flips. It should be noted that we do not take the 
global flips along the real-space direction into account in the present QMC calculation, assuming 
that their effect is not serious at least when the number N of spins is not too small. This means 
that we neglect the change of the winding number of the spin configuration.!^ 

Using this approach, we have investigated numerically the ground-state and thermodynamic 
properties of the spin-1/2-1/2-1-1 chain, confining ourselves to the case where J\ = J3 = 1.0 and 
J2 >0. (Note that Ji, or equivalently J3, is chosen to be the unit of energy.) 

3.2 Ground-state properties 

In this subsection, we discuss the energy gap and the nearest-neighbor two-spin correlation func- 
tions in the ground state, in order to investigate the ground-state phase transition. 

For a finite-iV system, we have carried out the QMC calculation only with the local flips. Thus, 
the system is kept in the subspace of a given value of the z-component 5 t 2 ot of <S to t- We nave 
estimated the energies E$* (N) with 5 t z ot = (the ground state) and 1 (the first-excited state) , 
which are expressed as 

E s * (N) = -(Qi) , (3.7) 

totV ' Y 6 /MC,S t z ot v ; 
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where 






N/4 \ ( p. (2r-l,2r) \ , /„(2r-l,2r) 

u PU-3.U-2 \ / (2r-1.2r) , / U P' ' 



r\ \ \ ./ I ' J M-3,U-2 I / (2r-l,2r) . / "HU-l.U \ / (2r-l,2r) 

Qi = 2^ U I — Q-p — / p «-3,«-2 + — ^ — I / pu-im 



r=l £=1 



(2r,2r+l) \ , /« (2r,2r+l)' 

(2r,2r+l) / °P4l,4l+l i , izr.zr-M, . ,_ * 

K4^-2,4£-l + f)0 I / r/«,4&fl f V°-°^ 



^41-2,41-1 \ / (2r,2r+l) / ^PuM+ 1 \ / (2r,2r+l) 

—QP / P4t-2,4i~l + 1 9(3 




and (• • ^mc.S 2 stands for the Monte Carlo average within the subspace determined by the value 
of 5 t z ot - The energy gap A(JV) is defined by 

A(N) = Ex(N) - E (N) . (3.9) 

Carrying out the above-mentioned QMC calculation within the 5 t z ot = subspace, we have also 
estimated the three kinds of ground-state nearest-neighbor two-spin correlation functions 

<*>i,2(N) = (si • s* 2 )mc,o = 3(sis|)mc,o , (3.10a) 

w 2 ,3(JV) = {82 ■ S 3 )mc,o = 3(s|S , |)mc,o , (3.10b) 

u; 3A (N) = (S 3 • S ? 4 )mc,o = 3(5|S|)mc,o • (3.10c) 

In these calculations, the values of J^ have been chosen to be J2= 1.00, 0.90, 0.80, 0.78, 0.75, 0.70, 
0.60, and 0.50, and those of N to be TV = 8, 16, 32, 64, and 128. We have performed 10 6 Monte Carlo 
steps after 10 5 initial steps for obtaining the thermal equilibrium. The Trotter numbers which we 
have used are n = 12, 16, 24, 32, 40, and 48, and the n-dependence of the finite- N QMC results 
A n (N) for a given physical quantity A has been extrapolated to n— >oo by making a least-squares 
fit to the formula^ 

A n (N)=A(N) + ^ + - i , (3.11) 

where a and b are constants which are independent of n; A(N) yields the Trotter-extrapolated 
value. Most of the calculations have been done at the temperature k-sT = 0.05 only. For the cases 
of J2 = 1.00 and 0.50 with N = 8 and 16, however, we have also done the calculation at k-gT = 0.02, 
and have obtained, for both temperatures, the same Trotter-extrapolated values for Eq(N) and 
Ei(N). These values agree with those obtained by the exact-diagonalization calculation within the 
numerical error, although the conservation of the winding number causes a systematic deviation in 
short chains. We therefore consider that k^T = 0.05 is low enough to discuss the zero-temperature 
properties. Thus, the calculated result for Eq(N) yields the ground-state energy E g (N) and those 
for u)i t 2(N), W2,3(iV), and uj 3t i(N) yield the nearest-neighbor two-spin correlation function in the 
ground state. 
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The Trotter-extrapolated values depend upon which data we use when we make a least-squares 
fit to the formula given by eq. (3.11). Following Miyashita and Yamamoto's procedure,™ we have 
performed three different extrapolations. The first one, denoted by TEX6, uses six data for n = 12, 
16, 24, 32, 40, and 48, the second one, TEX5, uses five data for n = 16, 24, 32, 40, and 48, and 
the third one, TEX4, uses four data for n = 24, 32, 40, and 48. Both for the ground-state energy 
eg(N) = Eo(N)/(N/4) per unit cell consisting of two 5 = 1/2 and two S = l spins and for the energy 
gap A(iV), the three extrapolations give almost the same values down to the second decimal place. 
We determine the Trotter-extrapolated value to be the average of the three values obtained by 
the three extrapolations, and estimate the error by the difference between the average and the 
farthest value among the three values. The results for e g (N) and A(iV) are tabulated, respectively, 
in Tables II and III, where, for the sake of comparison, the results of the exact-diagonalization 
calculation for N = 8 and 16 are also listed. 



Table II. Numerical results for the ground-state energy e s (N) per unit cell in the case where Ji = J3 = 1.0. The 
values in the columns denoted by TEX give the Totter-extrapolated results, and those in the columns denoted by 
ED give the results of the exact-diagonalization calculation. The extrapolated values e g (oo) to TV— +00 are given in 
the rightmost (TV— >oo) column. It is noted that the figures in the parentheses show errors in the last digit. 







N = 


= 8 


N = 


= 16 


7V = 32 


7V = 64 


TV = 128 


TV^oo 






TEX 


ED 


TEX 


ED 


TEX 


TEX 


TEX 




J 2 = 


= 1.00 


-3.54(1) 


-3.5470 


-3.52(1) 


-3.5159 


-3.52(1) 


-3.52(1) 


-3.52(1) 


-3.52(1) 


J 2 = 


= 0.90 


-3.37(1) 


-3.3832 


-3.34(1) 


-3.3449 


-3.34(1) 


-3.34(1) 


-3.34(1) 


-3.34(1) 


J 2 = 


= 0.80 


-3.22(1) 


-3.2326 


-3.19(1) 


-3.1891 


-3.18(1) 


-3.18(1) 


-3.18(1) 


-3.18(1) 


J 2 = 


= 0.78 


-3.19(1) 


-3.2045 


-3.16(1) 


-3.1609 


-3.15(1) 


-3.15(1) 


-3.15(1) 


-3.15(1) 


J 2 = 


= 0.75 


-3.15(1) 


-3.1638 


-3.12(1) 


-3.1209 


-3.11(1) 


-3.11(1) 


-3.11(1) 


-3.11(1) 


J 2 = 


= 0.70 


-3.09(1) 


-3.1002 


-3.06(1) 


-3.0611 


-3.05(1) 


-3.05(1) 


-3.05(1) 


-3.05(1) 


J 2 = 


= 0.60 


-2.98(1) 


-2.9908 


-2.97(1) 


-2.9655 


-2.96(1) 


-2.96(1) 


-2.96(1) 


-2.96(1) 


J 2 = 


= 0.50 


-2.90(1) 


-2.9060 


-2.89(1) 


-2.8937 


-2.89(1) 


-2.89(1) 


-2.89(1) 


-2.89(1) 



In order to obtain the results in the thermodynamic limit (iV — >oo), the Trotter-extrapolated 
values A(N) have further been extrapolated by the least-squares method using a linear function of 
l/N 2 , 

A(N) = A(oo) + -^ (3-12) 
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Table III. Numerical results for the energy gap A(iV) in the case where Ji = J3 = 1.0. The values in the columns 
denoted by TEX give the Totter-extrapolated results, and those in the columns denoted by ED give the results of 
the exact-diagonalization calculation. The extrapolated values A(oo) to N—>oo are given in the rightmost (N— +00) 
column. It is noted that the figures in the parentheses show errors in the last digit. 







N = 


= 8 


N = 


= 16 


iV = 32 


iV = 64 


iV = 128 


JV-fOo 






TEX 


ED 


TEX 


ED 


TEX 


TEX 


TEX 




J2-- 


= 1.00 


0.61(1) 


0.6156 


0.46(1) 


0.4497 


0.44(1) 


0.44(1) 


0.44(1) 


0.44(1) 


J2-- 


= 0.90 


0.55(1) 


0.5566 


0.36(1) 


0.3522 


0.32(1) 


0.31(1) 


0.31(1) 


0.31(1) 


J2-- 


= 0.80 


0.51(1) 


0.5126 


0.28(1) 


0.2807 


0.17(1) 


0.13(1) 


0.13(1) 


0.13(1) 


J2-- 


= 0.78 


0.50(1) 


0.5064 


0.27(1) 


0.2733 


0.15(1) 


0.09(1) 


0.06(2) 


0.06(2) 


J2-- 


= 0.75 


0.50(1) 


0.4992 


0.27(1) 


0.2685 


0.14(1) 


0.08(1) 


0.06(2) 


0.06(2) 


J2-- 


= 0.70 


0.49(1) 


0.4930 


0.28(1) 


0.2785 


0.20(1) 


0.18(1) 


0.18(1) 


0.18(1) 


J2-- 


= 0.60 


0.51(1) 


0.5059 


0.37(1) 


0.3555 


0.34(1) 


0.34(1) 


0.34(1) 


0.34(1) 


J2-- 


= 0.50 


0.56(1) 


0.5529 


0.48(1) 


0.4697 


0.48(1) 


0.48(1) 


0.48(1) 


0.48(1) 



with a constant c. Performing the extrapolation, we have used three data 4(32), 4(64), and 4(128) 
for J 2 = 0.75 and 0.78, and four data 4(16), 4(32), 4(64), and 4(128) for the other values of J 2 . 
Figure 4 illustrate this extrapolation for the energy gap, where 4 is A. The N^oo extrapolated 
values e g (oo) for the ground-state energy per unit cell and those A (00) for the energy gap are also 
listed in Tables II and III, respectively, and are plotted as a function of J 2 in Fig. 5 and Fig. 6, 
respectively. From Fig. 6 we see that A(oo) vanishes at 



J-> = J: 



2c 



0.77 ±0.01 



(3.13) 
1.0 undergoes a 



This result shows that the ground state of the present system with J\ = J3 
second-order phase transition at this value of J 2 . 

The Trotter extrapolated values ur i2 (iV), w 2i 3(-/V), and 0^4 (iV) for the nearest-neighbor two-spin 
correlation functions are tabulated in Tables IV, V, and VI, respectively, where the corresponding 
N — > ex) extrapolated values, wi i2 (oo), u; 2j 3(oo), and 0^4(00) are also listed. In Fig. 7 we plot 
tJi j2 (oo), w 2i 3(oo), and ^3,4(00) as functions of J 2 , where each correlation function is normalized 
by its minimum value (note that ~3/4<wi )2 (iV)<l/4, -l<W2,3(JV)<l/2, and -2<w 3i4 (iV) < lY 
This figure demonstrates that as J 2 increases, |u;i !2 (oo)| and 1^2,3 (00) | rapidly decreases and in- 
creases, respectively, around J2~J2c, while |o;3 i 4(oo)| rather gradually decreases with increasing 
J 2 . Giving a thought to these J 2 -dependences of the correlation functions, we may schematically 
represent, by means of the VBS picture,^ the ground states for J 2 < J 2c and that for J 2 > J 2c 
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as depicted in Fig. 8. The VBS wave function of Fig. 8(a) corresponds to the state for J 2 = 0.0, 
while the VBS wave function of Fig. 8(b) to the state for Ji = 0.0 and J2^+oo. The change of 
the correlation functions as a function of J2 clearly illustrates the interpolation between these two 
limiting states. 



Table IV. Trotter-extrapolated ground-state nearest-neighbor two-spin correlation function u>i.2(N) and the TV^oo 
extrapolated ground-state nearest-neighbor two-spin correlation function 011,2(00) in the case where Ji = J3 = 1.0. 
It is noted that the figures in the parentheses show errors in the last digit. 





iV = 8 


iV = i6 


iV = 32 


JV = 64 


iV = 128 


TV^oo 


J 2 = 1.00 


-0.28(1) 


-0.23(1) 


-0.22(1) 


-0.22(1) 


-0.22(1) 


-0.22(1) 


J 2 = 0.90 


-0.33(1) 


-0.27(1) 


-0.26(1) 


-0.26(1) 


-0.26(1) 


-0.26(1) 


J 2 = 0.80 


-0.41(1) 


-0.37(1) 


-0.35(1) 


-0.34(1) 


-0.34(1) 


-0.34(1) 


J 2 = 0.78 


-0.43(1) 


-0.40(1) 


-0.39(1) 


-0.38(1) 


-0.38(1) 


-0.38(1) 


J 2 = 0.75 


-0.46(1) 


-0.46(1) 


-0.46(1) 


-0.46(1) 


-0.46(1) 


-0.46(1) 


J 2 = 0.70 


-0.51(1) 


-0.53(1) 


-0.54(1) 


-0.54(1) 


-0.54(1) 


-0.54(1) 


J 2 = 0.60 


-0.60(1) 


-0.62(1) 


-0.62(1) 


-0.62(1) 


-0.62(1) 


-0.62(1) 


J 2 = 0.50 


-0.66(1) 


-0.67(1) 


-0.67(1) 


-0.67(1) 


-0.67(1) 


-0.67(1) 



3.3 Thermodynamic properties 

We now turn to the discussion of the thermodynamic properties. Performing the QMC calculation 
which takes into account the global flips given by eqs. (3.5a) and (3.5b) as well as the local flips 
given by eqs. (3.5a)-(3.5e), we have calculated the temperature dependences of the specific heat 
C(N) and the magnetic susceptibility x(iV) for J 2 = J2 C = 0.77, and also for J2 = 1.0 and 5.0. Here, 
the expressions for C(N)£y and x(iV) are given by 



C(N) 




MC 



Ch 



and 



X(N) 




MC \ /MC 



MC 



(3.14) 



(3.15) 



where (• • -)mc denotes the Monte Carlo average, Qi is given by eq. (3.8), and 
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Table V. Trotter-extrapolated ground-state nearest-neighbor two-spin correlation function u)2,3(N) and the N^oo 
extrapolated ground-state nearest-neighbor two-spin correlation function 0^2,3(00) in the case where Ji = J3 = 1.0. 
It is noted that the figures in the parentheses show errors in the last digit. 





iV = 8 


JV = 16 


iV = 32 


JV = 64 


iV = 128 


N^oo 


J 2 = 1.00 


-0.87(1) 


-0.89(1) 


-0.89(1) 


-0.89(1) 


-0.89(1) 


-0.89(1) 


J 2 = 0.90 


-0.81(1) 


-0.84(1) 


-0.85(1) 


-0.85(1) 


-0.85(1) 


-0.85(1) 


J 2 = 0.80 


-0.72(1) 


-0.74(1) 


-0.75(1) 


-0.76(1) 


-0.76(1) 


-0.76(1) 


J 2 = 0.78 


-0.70(1) 


-0.70(1) 


-0.71(1) 


-0.72(1) 


-0.72(1) 


-0.72(1) 


J 2 = 0.75 


-0.67(1) 


-0.64(1) 


-0.63(1) 


-0.62(1) 


-0.62(1) 


-0.62(1) 


J 2 = 0.70 


-0.61(1) 


-0.54(1) 


-0.52(1) 


-0.51(1) 


-0.51(1) 


-0.51(1) 


J 2 = 0.60 


-0.47(1) 


-0.41(1) 


-0.40(1) 


-0.40(1) 


-0.40(1) 


-0.40(1) 


J 2 = 0.50 


-0.35(1) 


-0.31(1) 


-0.31(1) 


-0.31(1) 


-0.31(1) 


-0.31(1) 



Table VI. Trotter-extrapolated ground-state nearest-neighbor two-spin correlation function u>3,4(iV) and the N -^00 
extrapolated ground-state nearest-neighbor two-spin correlation function 013,4(00) in the case where Ji = J3 = 1.0. 
It is noted that the figures in the parentheses show errors in the last digit. 





iV = 8 


iV = 16 


iV = 32 


JV = 64 


iV = 128 


N^oo 


J 2 = 1.00 


-1.58(1) 


-1.53(1) 


-1.52(1) 


-1.52(1) 


-1.52(1) 


-1.52(1) 


J 2 = 0.90 


-1.62(1) 


-1.57(1) 


-1.56(1) 


-1.56(1) 


-1.56(1) 


-1.56(1) 


J 2 = 0.80 


-1.69(1) 


-1.66(1) 


-1.64(1) 


-1.63(1) 


-1.63(1) 


-1.63(1) 


J 2 = 0.78 


-1.70(1) 


-1.68(1) 


-1.67(1) 


-1.66(1) 


-1.66(1) 


-1.66(1) 


J 2 = 0.75 


-1.73(1) 


-1.73(1) 


-1.72(1) 


-1.72(1) 


-1.72(1) 


-1.72(1) 


J 2 = 0.70 


-1.78(1) 


-1.79(1) 


-1.80(1) 


-1.80(1) 


-1.80(1) 


-1.80(1) 


J 2 = 0.60 


-1.86(1) 


-1.87(1) 


-1.87(1) 


-1.87(1) 


-1.87(1) 


-1.87(1) 


J 2 = 0.50 


-1.92(1) 


-1.92(1) 


-1.92(1) 


-1.92(1) 


-1.92(1) 


-1.92(1) 



Qi 



n N/4 

EE 

r=l 1=1 



pa (2r-l,2r) 
u P4^-3,4^-2 



(2r- l,2r) 
^4£-3,4^-2 



p. (2r-l,2r) 

dp 



(2r-l,2r) 
PU-3,U-2 
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8 2 (2r-l,2r)\ , /- (2r-l,2r) \ / 

^ F4l-l,4l 1 / (2r-l,2r) _ I / ^Alg-l^l | / (2r-l,2r) 



r,2 (2r,2r+l) \ / | /o (2r,2r+l) 

^ ^41-2,41-1 \ / (2r,2r+l) _ I [ a PU-2,U-l \ I (2r,2r+l) 
+ I Qfi2 I / ^4£-2,4£-l ~ S I go I / PU-2M-1 



8 2(2r,2r+l)\ , f / „ (2r,2r+l) \ , 

/ a K4l,4l+1 I / (2r,2r+l) _ I / a F4l,4l+l \ / (2r,2r+l) 
+ ! 9/3 2 I / P«,«+l "SI ^ I / ^,4£+l 



(3.16) 



In deriving eq. (3.15), we assume that the g-factor associated with the 15* = 1/2 spin is identical to 
that associated with the S = l spin and that (5^ ot )MC vanishes. For J2 = 1.0 and 5.0, the values of 
iV have been chosen only to be N = 8, 16 and 32. This is because for these J2's, we expect weak 
iV-dependences of C(N)/(N/4) and x(iV)/(JV/4), since the energy gap A(oo) is rather large (see 
Fig. 6), or equivalently, the correlation length is rather short. For the critical value J2 = 0.77, on 
the other hand, the calculation has been done also for N = 64 and 128 when T is not too high, in 
addition to N = 8, 16 and 32. We list in Table VII Trotter numbers n's, the Monte Carlo steps, and 
the initial steps for obtaining the thermal equilibrium. By using the QMC results for all of the four 
Trotter numbers, the n-dependence of C n {N) for J2 = 0.77 and 1.0 has been extrapolated to ra— >oo 
by the least-squares method with the formula of eq. (3.11)™) Similar extrapolations for C n (N) for 
J2 = 5.0 and for Xn(N) for J2 = 0.77, 1.0, and 5.0 have been done by the use of the formula© 

A n (N) = A(N) + ^ , (3-17) 

where a is a constant, instead of the formula of eq. (3.11). 

In Fig. 9 we plot versus the temperature T the Trotter-extrapolated values C{N)/(N/4) [Fig. 9(a)] 
with N = 8, 16, and 32 for the specific heat per unit cell and x(A r )/(A r /4) [Fig. 9(b)] with N = 8, 
16, and 32 for the magnetic susceptibility per unit cell, which have been obtained for J2 = 1.0. For 
the sake of comparison, we also plot, in these figures, the exact results for both quantities for iV = 8 
and J2 = 1.0 versus T, which we have obtained by diagonalizing completely the Hamiltonian to 
calculate all the eigenvalues. The corresponding plots for J2 = 0.77 are presented in Fig. 10. As is 
expected, both C{N)/{N/4) and x(-^)/(-^/4) for J2 = 1.0 have rather weak iV-dependences even 
at /cBr = 0.08. For J2 = 0.77, on the other hand, the iV-dependence of x{N)/{N/A) is noticeable 
at low temperatures; it increases as N increases. Judging from the T-dependences of x(64)/16 
and x(128)/32, we may conclude that at T = 0.0, the magnetic susceptibility per unit cell in the 
thermodynamic limit is finite for J2 = 0.77 with the value [x(A'")/(A^/4)]jv_^ 00 =0.42 ± 0.01. This is 
consistent with the fact that the energy gap vanishes for this value of Ji- 
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Table VII. Sets of the Trotter numbers (n) used at various temperatures, and also the Monte Carlo steps (MCS) 
and the initial steps (IS) spent for each Trotter number at each temperature, where the former steps do not include 
the latter steps. Note that J1—J3 — I.O. 

k B T for J 2 = 0.77, 1.0 fc B TforJ 2 = 5.0 n MCS IS 



0.08<fc B T<0.10 


0.400 <k B T< 0.625 


12, 16, 24, 32 


120 xlO 4 


6 xlO 4 


fc B T = 0.15 


fc B T = 0.750 


8, 12, 16, 24 


100 xlO 4 


5 xlO 4 




fc B T = 0.875 


8, 12, 16, 24 


50xl0 4 


5 xlO 4 


fc B T = 0.20 


1.000 <k B T< 1.375 


6, 8, 12, 16 


50 xlO 4 


5 xlO 4 


0.30 <k B T< 0.50 


1.500 <k B T< 2.750 


4, 6, 8, 12 


40 x 10 4 


5 xlO 4 


0.60 <k B T< 1.00 


3.000 <k B T< 5.000 


2, 4, 6, 8 


30 xlO 4 


5 xlO 4 


1.20 <k B T< 2.00 


6.000 <k B T< 8.000 


2, 4, 6, 8 


25 xlO 4 


5 xlO 4 



Figure 11 presents the plots versus T of C{N)/(N/A) [Fig. 11(a)] with 7V = 8, 16, and 32 and 
x(iV)/(iV/4) [Fig. 11(b)] with N = 8, 16, and 32, which have been obtained for J2 = 5.0, and also 
the plots versus T of the exact results for both quantities for N = 8 and J2 = 5.0. It is noted that 
the iV-dependences of both C{N)/{N/A) and x(iV)/(iV/4) are weak as in the case of J2 = 1.0. A 
characteristic feature of the C{N)/{N/A) versus T curve shown in Fig. 11(a) is that it has a double 
peak. The origin of the two peaks can be understood as follows. The dot-dashed line in Fig. 11(a) 
shows the T-dependence of the specific heat for two pairs of the two-spin system consisting of an 
S=l/2 spin s and anS=l spin S which couple with each other by 5s ■ S. Comparing this with the 
T-dependence of C(iV)/(iV/4), we may conclude that the higher-temperature peak is associated 
with the Schottky-type peak of this two-spin system. On the other hand, the dashed line in 
Fig. 11(a) shows the T-dependence of the half of the specific heat for the S= 1/2 bond-alternating 
chain described by the Hamiltonian T~L< a ) [ see e Q- (2-3)] with J\ = J3 = 1.0 and N = 8. The peak 
height as well as the peak position in this T-dependence agrees, respectively, with those in the 
T-dependence of C(N)/(N/4), which may leads us to the conclusion that the lower-temperature 
peak is associated with the peak coming from the short-range order in the 5=1/2 bond-alternating 
chain. A somewhat unusual behavior of x(7V)/(iV/4) as a function of T, seen in Fig. 11(b), can also 
be understood as follows. At high temperatures the T-dependence of x(iV)/(iV/4) should be similar 
to that of the magnetic susceptibility for two pairs of the above two-spin system, which diverges 
at T^O.O [see the dot-dashed line in Fig. 11(b)], while x(iV)/(iV/4) should vanish at T = 0.0 since 
the present system has massive excitations when J2 = 5.0. 



18 



§4. Summary and Conclusions 

We have explored the ground-state properties of the quantum mixed spin system described by the 
Hamiltonian TC see eq. (1.1) • First, we have applied the Lieb-Mattis theoremliZP to show that the 
ground state of the system is nonmagnetic when J% > and when J\ < and J3 > 0; in the remaining 
region where J\ < and J3 < it is ferrimagnetic or ferromagnetic depending upon whether J2 > 
or J2<0 Table I . Second, applying the Wigner-Eckart theoremO- 1 to the limiting cases (a) 
I J2 1 ^ \Ji I > I <^3 1 with J2 > 0, (b) I J2 1 S> I J\ I , I J3 1 with J2 < 0, and (c) | Ji | 3> | J2 1 , \J$\ with J\ < 0, we 
have mapped the Hamiltonian 7i into the bond-alternating chains described by the Hamiltonians 



Ti.( a ) with effective 5 = 1/2 spins 
with effective S = 1 spins eq. (2.7) 



eq. (2.3) , ~H{p) with effective S = 3/2 spins eq. (2.5) , and T~i{ c ) 
, respectively. This leads to the facts that in the limiting case (a) 
the present system is massless when eq. (2.8) holds, and in the limiting case (b) it is massless when 
eq. (2.9) holds. We have also carried out the perturbation calculations to show that in the limiting 
cases (d) J3» I Ji|, \Ji\ and (e) J\ 3> | J2I and J\ ^> J3 >0, the system is massless when eqs. (2.14) and 
(2.18) hold, respectively. These results for J3 = 1.0 are summarized in Fig. 2. Third, performing the 
QMC calculation without the global flips at a sufficiently low temperature (A;bT' = 0.05), we have 
shown that when Ji = J^ = 1.0, the energy gap A(oo) vanishes at J2 = J2 C = 0.77 ± 0.01 Fig. 6 , 
where the second-order phase transition occurs in the ground state. The representation of both 
the ground state for J2 < J2 C and that for J2 > J<i c by means of the VBS picture is given in Fig. 8, 
which is suggested from the J2 dependences Fig. 7 of the ground-state nearest-neighbor two-spin 
correlation functions 11^2(00), 1^2,3(00), and 0^4(00). 

Furthermore, performing the QMC calculation including the global flips along the Trotter direc- 
tion as well as the local flips, we have calculated the temperature dependences of the specific heat 
C{N)/{N/4) per unit cell and the magnetic susceptibility x(iV)/(iV/4) per unit cell with N = 8, 
16, and 32 for J 2 = 0.77, 1.0, and 5.0 with J\ and J3 fixed at Ji = J% = 1.0. The results of these 
calculations are depicted in Figs. 9-11. 

Figure 2 clearly demonstrates that several massless lines exist on the J\ versus J2 plane with 
J\ >0 in the case of J3 = 1.0. We are now exploring the details of the ground-state phase diagram 
on this plane by the use of a density-matrix renormalization-group method proposed originally 
by WhiteBa) Our preliminary results show that there exist four massless lines which divides the 
upper-half plane into six regions; the ground state in each region can be understood by the VBS 
picture.™ 

In conclusion, we hope that the present study stimulates further experimental studies on related 
subject, which include the synthesization of quantum mixed spin systems with nonmagnetic ground 
states. 
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Fig. 1. Unit cells of the lowest-energy spin arrays, for (a) J3 >0 and (b) J3 <0, of the corresponding classical system 
in which at and St are replaced, respectively, by the classical spin vectors with the magnitudes of 1/2 and 1. The 
short and long arrows indicate the former and latter classical spin vectors, respectively. 



Fig. 2. Massless lines in the J\ versus J2 plane with Ji >0 and J3 = 1.0 in the limiting cases (a), (b), (d), and (e). 
The massless lines in the cases (a), (b), and (e) are shown by the dashes lines, while the massless line in the case 
(d) is by the solid line in the box where the scales of both the abscissa and the ordinate are enlarged by factor 
10. The dotted lines show simple extrapolations of the solid line. The cirle denotes the massless point (J2 =0.77, 
Ji =1.0), which is obtained in §3.2 [see eq. (3.13)]. 



Fig. 3. A graphical representation of the two-dimensional Ising system, where the horizontal and vertical directions 
correspond, respectively, the real-space and Trotter directions. The open circles denote the Ising variables {s^_ 3 } 
and {s^_ 2 }, which take the two values ±^, and the open squares denote the Ising variables {S^,^}, and {S^ }, 
which take the three values and ±1. The plaquettes shaded by the slashed, back-slashed, vertical, and horizontal 
lines correspond, respectively, to the local Boltzmann factors, p 4 [Sz ie-2> Pkt-i kt > Pk[-2M-n an d Pktkt+i ■ 



Fig. 4. Plot of the Trotter-extrapolated energy gap A(N) versus 1/N with N= 16, 32, 64, and 128, for (a) J 2 = 1.00 
(the open circles), J2 = 0.90 (the open squares), J2=0.80 (the open diamonds), and J2=0.78 (the open triangles) 
and for (b) J2 =0.75 (the closed triangles), J2 =0.70 (the closed diamonds), J2 =0.60 (the closed squares), and 
J2=0.50 (the closed circles), where Ji = J3 = 1.0. The solid lines give curves of a least-squares fit to eq. (3.12) (see 
the text for more details). 
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Fig. 5. Plot versus J2 of the N —* 00 extrapolated ground-state energy per unit cell e g (oo) , where 3\ = J3 = 1.0. 
The solid line is a guide to the eye. 



Fig. 6. Plot versus J2 of the N — > 00 extrapolated energy gap A(oo), where J\ = J3 = 1.0. The solid lines are guides 
to the eye. 



Fig. 7. Plots versus J2 of the iV— >co extrapolated ground-state nearest-neighbor two-spin correlation functions 
0*1,2(00) (the open circles), 012,3(00) (the open squares), and 0)3,4(00) (the open diamonds), each of which is nor- 
malized by its minimum value. Note that J\ — J3 = 1.0. The solid lines are guides to the eye. 



Fig. 8. Schematic representations of the ground states for (a) J2 < Jic and (b) J2 > J2 C by means of the VBS 
picture£j' where J1 — J3 — 1.0. The solid circles represent the S— 1/2 spins, and two S — 1/2 spins connected by 
the solid line form a singlet pair. Each open ellipse surrounding two 5=1/2 represents an operation of constructing 
an S=l spin from these S— 1/2 spins by symmetrizing them. 



Fig. 9. Plot versus the temperature T of (a) the Trotter-extrapolated specific heat C(N)/(N/4) per unit cell and 
(b) the Trotter-extrapolated magnetic susceptibility ^(JV)/(JV/4) per unit cell with N = 8 (the open circles), 16 
(the crosses), and 32 (the open diamonds) for J2 = 1.0. The dotted lines in (a) and (b) show, respectively, the exact 
results for the specific heat per unit cell and the magnetic susceptibility per unit cell, for N = 8 and J2 = 1.0. Note 
that Ji = J 3 = 1.0. 
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Fig. 10. Plot versus the temperature T of (a) the Trotter-extrapolated specific heat C(N)/(N/4) per unit cell and 
(b) the Trotter-extrapolated magnetic susceptibility x(N)/(N/4:) per unit cell with N — 8 (the open circles), 16 (the 
crosses), and 32 (the open diamonds), 64 (the open squares), and 128 (the open triangles) for J2=0.77, being the 
critical value of J2. Here, C(N) and x(N) with iV = 64 and 128 are given only when fcBT<0.50. The dotted lines 
in (a) and (b) show, respectively, the exact results for the specific heat per unit cell and the magnetic susceptibility 
per unit cell, for iV = 8 and J 2 = 0.77. Note that Ji = J 3 = 1.0. 



Fig. 11. Plot versus the temperature T of (a) the Trotter-extrapolated specific heat C(N)/(N/4) per unit cell and 
(b) the Trotter-extrapolated magnetic susceptibility x{N)/(N/&) per unit cell with N — 8 (the open circles), 16 
(the crosses), and 32 (the open diamonds) for J2 =5.0. The dotted lines in (a) and (b) show, respectively, the 
exact results for the specific heat per unit cell and the magnetic susceptibility per unit cell, for N = 8 and J2 = 5.0. 
Note that J\ — Js — 1.0. The dot-dashed lines in (a) and (b) show, respectively, the specific heat and the magnetic 
susceptibility for two pairs of the two-spin system consisting of an S=l/2 spin s and an5 = l spin S which couple 
with each other by 5s • S. Furthermore, the dashed line in (a) shows the half of the specific heat for the S— 1/2 
bond-alternating chain described by the Hamiltonian 7~C( a ) given by eq. (2.3) with Ji — J 3 — 1.0 and N = 8. 
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